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The Lya forest has emerged as one of the few systems capable of probing small-scale structure 
at high-z with high precision. In this talk I highlight two areas in which the Lya forest is shedding 
light on fundamental questions in cosmology, one speculative and one which should be possible in 
the near future. 
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I. INTRODUCTION 

As we heard in the talks by Uros Seljak and Lam Hui 
in this meeting, the physics of the Lya forest, as mea- 
sured in high-z Quasi-Stellar Object (QSO) spectra, is 
reasonably simple. At these redshifts the gas making up 
the intergalactic medium is in photoionization equilib- 
rium, which results in a tight density-temperature rela- 
tion for the absorbing material with the neutral hydrogen 
density proportional to a power of the baryon density. 
Since pressure forces are sub-dominant, the neutral hy- 
drogen density closely traces the total matter density on 
the scales relevant to the forest (0.1 — 10 7i _1 Mpc). The 
structure in QSO absorption thus traces, in a calcula- 
ble way, slight fluctuations in the matter density of the 
universe back along the line of sight to the QSO, with 
most of the Lya forest arising from over-densities of a 
few times the mean density. 

In this contribution I want to focus on two areas where 
the Lya forest can be used to constrain cosmology and 
astrophysics. The first, dealing with measuring the mat- 
ter power spectrum at high redshift, is quite speculative, 
but might be very exciting. The second has a higher 
chance of success and may tell us something about the 
nature of the ionizing sources, the end of the dark ages 
and the formation of the first objects in the universe. 
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FIG. 1: The contribution to the dark matter variance per 
logfc, A 2 (fc), in linear theory. The upper panel shows A 2 (fc) 
on a logarithmic scale while the lower panel shows the ratio 
of the numerically computed A 2 (fc) to the fitting function 
of |Tl| with the oscillations suppressed. 



II. FINDING BARYONS AT z > 1? 

The presence of a series of almost regularly spaced 
peaks in the cosmic microwave background tempera- 
ture angular power spectrum has now been compellingly 
demonstrated. The same physics which gives rise to these 
peaks, acoustic oscillations in the baryon-photon plasma 
prior to recombination , also predicts a (much smaller) 
series of peaks in the matter power spectrum. A mea- 
surement of these peaks would provide a fundamental 
test of the paradigm, and an independent measure of the 
z-dependent Hubble constant and distance-redshift rela- 
tion 0, 0, 0, IE IE ■ Because non-linear growth 'washes 
out' the peaks |E they are best searched for at high red- 
shift. 

The baryon oscillations themselves lie at k ~ 0.03 — 
0.3/i/Mpc (see Fig. QJ, with an amplitude of 5% or less 
for currently popular cosmological models 0, 0|. The 
higher modes are roughly out of phase with the pho- 



ton peaks (also in fc-space) and appear half as often 
as the peaks in the CMB spectrum (see e.g. Fig. 1 of 
0). The peaks are suppressed by Silk damping beyond 
k ~ 1 h/Mpc. 

Naively the best way to measure the small baryonic 
oscillations in the matter power spectrum would be an 
ultra-large galaxy redshift survey at high-z. Such a 
survey would require next-generation instruments and a 
large dedicated team. Thus even with adequate funding 
and heroic efforts it is several years away. It is there- 
fore interesting to explore whether one could measure 
the baryonic oscillations another way. 

Within a few years the SDSS will return a catalogue 
of ~ 10 5 QSOs over ~ 10 4 square degrees of sky. These 
quasars will be distributed in both redshift and angular 
separation, but simple counting suggests that one should 
find many hundreds to several thousands with separa- 
tions of tens of Mpc above z ~ 2. It is therefore inter- 
esting to ask whether this sample could be used to probe 
the baryon oscillations. 
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Using simply the auto-power spectrum of the forest 
will not work, because the line-of-sight spectrum at k is 
an integral of the 3D spectrum over wavenumbers greater 
than k. This integration suppresses the small oscillations 
which are our primary signal. However one can con- 
sider the cross-spectra, where one line-of-sight is cross- 
correlated with a neighboring line-of-sight. As we shall 
show below this statistic is sensitive to the oscillations. 
In addition uncorrelated errors from e.g. continuum fit- 
ting, should be averaged down in this procedure relative 
to their value in the auto-spectrum. Since we are work- 
ing on large scales, where the noise is essentially absent, 
one could in principle measure the flux cross-spectrum to 
the cosmic variance limit 



5n T ^ / 2k f 
77 ~ V NAk 



(1) 



where N is the number of pairs, kf is the fundamental 
mode along the line-of-sight and Ak is the width of the 
bin in fc-space in which we are measuring the power. We 
can alternatively write kt = 2tt/L, where L is the length 
of the 'useful' part of the spectrum. Unless we attempt to 
fit to the higher lines, the useful length of the spectrum is 
the section between Lya and Ly/3 which covers hundreds 
of Mpc in the rest frame of the absorbers. Thus with 
a few hundred pairs it is possible to measure the cross- 
spectrum with error bars on the order of per cent. 

One major advantage here is that, like the CMB itself, 
the redshift and scale dependence of the peaks is accu- 
rately calculable given a cosmological model. This allows 
us to perform a likelihood estimation of the relevant cos- 
mological factors, including pairs of different distances 
and redshifts, rather than attempt to perform an inver- 
sion on the irregularly sampled data. Though the latter 
can be done, it is considerably more difficult than the 
'forward' approach. 

By applying the idea of the peak-background split from 
large-scale structure is it possible to show that, for suffi- 
ciently large scales and assuming r = A(p/p) a , the flux 
cross-spectrum, 7r r , along two lines of sight separated 
by distance x± is proportional to the dark matter cross- 
power spectrum it (see also On these very large 
scales we can neglect the effects of thermal broadening 
or peculiar velocities. The DM cross-spectrum is simply 
related to the 3D power spectrum by 



TT 2 (k,X ± ) = 



n d -i ( k -\ u^k2)A\ q ) ( 2 ) 

Jk <7 \1/ 



or 



m 



A 2 (<z) = / 



dx 



(qx) 2 Jo(xVq^¥)TT 2 (q) (3) 



where the inverse relation holds V k < q. Note that for 
smooth spectra, 7r 2 (fc) gets contributions from k < q < 
xj_ . As x± — > 0, Jo — > 1 and we find the auto-spectrum 
is the projection of the 3D spectrum as anticipated. 
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FIG. 2: The cross-power spectra of the DM, in linear theory, 
at z = 3. The upper panel shows n(k,x±), including and 
excluding the baryonic oscillations, for different values of x±. 
The lower panel shows the ratio to the smooth fitting function 

of EH. 



Converting from h/Mpc to s/km using H(z) = 
H yJVL M (l + zf + Qa ~ lOOkm/s/Mpc we find that 
the relevant spectral range is roughly 3 x 10~ 4 s/km to 
3 x 10~ 3 s/km. This is quite large scale, roughly compa- 
rable to the fundamental mode if we measure only from 
lOOnm to 120nm. 

On the positive side Fig. [2] shows that the signal is 
present in the cross spectra at separations of lOMpc or 
so. The relative amplitude of the oscillations (which is 
the figure of merit for cosmic variance dominated mea- 
surements) can be quite large, 10% or more, compared 
to the accuracy with which it can be measured. Because 
we can accurately predict the positions and amplitudes 
of the peaks we know both the x± and z dependence of 
the signal. 

On the negative side there is only a small fractional 
change in ir when it is big (i.e. x± small) and the larger 
fractional change in ir when ir is small. This suggests 
we would be working primarily at very low signal levels, 
where systematics become increasingly important. Also, 
if we work with the spectra only between Lya and Ly/3 
we do not have very much 'spare' resolution in k to see 
the oscillations. We would ideally probe k <C 0.003s/km 
which requires very long spectra. 

In summary, the signal of baryon oscillations is lurk- 
ing in the cross-spectra of QSO pairs, which allow us 
to probe large volumes of space at high redshift with a 
modest investment in telescope time. A naive estimate of 
the cosmic variance limitations of a QSO survey like the 
SDSS suggests that the signal could be seen in a likeli- 
hood analysis with high confidence. However the overall 
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signal levels are very low and there are numerous sources 
of systematic error which may make the signal forever 
unmeasurable. 



III. THE NATURE OF THE IONIZING 
SOURCES 



Currently the sources that make up the UV back- 
ground are unknown. However the mean optical depth 
f in the Lya forest offers a way to constrain the source 
emissivity and fluctuations in the UVbg can be used to 
constrain the nature of the sources. 1 have been involved 
in a long-running collaboration with Avery Meiksin to 
investigate these issues 0, 113 , and present some of our 
basic results here. 

First we consider the constraints on the source emis- 
sivity e. Given a cosmological model, the mean flux or 
optical depth f fixes the ionization rate T. If we approx- 
imate exp[— r] as exp[— £/ro] where ro is the attenuation 
length then T ~ er c g. Here r ff is an effective attenua- 
tion length which interpolates between the attenuation 
limited (r <C In) and cosmological expansion+age lim- 
ited (r ^S> £h) regimes and £r is the Hubble distance. 

As a function of e the transition between the attenu- 
ation limited and cosmologically limited regimes is very 
abrupt, and at fixed e, ro falls rapidly with increasing 
redshift. This in turn means that f rises rapidly, almost 
exponentially, with redshift. 

Fig. shows the emissivity of the sources, as a func- 
tion of z, required to fit the observed f. This is compared 
to extrapolations of the QSO luminosity function. From 
this we can see that QSOs may contribute a substantial 
fraction of the UV background even as high as z ~ 6. 
If QSOs do contribute non-negligibly to the UV back- 
ground then we would expect large spatial fluctuations 
in r. Conversely, if galaxies or star clusters dominate 
the UV background we expect to see only the imprint 
of large-scale structure. Unfortunately the effects of UV 
background correlations are most conspicuous at z > 5, 
where the precise level of absorption by the Lya forest is 
most difficult to measure because of extremely low flux 
values. 

Fig.0|shows how a simulated Lya spectrum is changed 
if the UV background is (unrealistically) dominated by a 
single source as compared to the usually assumed smooth 
behavior. To assess a more realistic situation it is nec- 
essary to resort to Monte Carlo simulations in which 
sources are placed within the simulation volume, the UV 
background computed and then spectra created from the 
density and velocity fields within the box. Our inves- 
tigations suggest that the 1-point distribution is rela- 
tively unaffected by fluctuations in the UVbg, however 
the 2-point statistics (the flux correlation function and 
power spectrum) show noticeable differences at high red- 
shift [LJ . In addition fluctuations in the UVbg will in- 
crease the estimated Lya optical depth over the value 
that would be estimated assuming a homogeneous back- 
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FIG. 3: (Upper) The emissivity, at the Lyman edge, required 
to fit the mean flux observed in QSO spectra as a function 
of wavelength (points) compared to extrapolated QSO emis- 
sivities for 3 different assumptions about the evolution of the 
QSO luminosity function (lines). Under these assumptions 
QSOs would be enough to keep the universe ionized at z < 6. 
(Lower) The attenuation length, ro in proper Mpc/ft, as a 
function of z. Note that ro drops quasi-exponentially to z ~ 6. 
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FIG. 4: (Upper) The normalized flux, as a function of veloc- 
ity, from a PM simulation of structure formation. The spectra 
were computed as described in [T3,ll3j assuming either a uni- 
form ionizing background (dashed) or that the background 
was dominated by a single source at the center of the box 
(dotted). (Lower) The dark matter density (dotted) along 
the line-of-sight or the inverse ionization rate, T -1 , in units 
of the average (dashed). 
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ground |12(. This increases the demand placed on the 
ionization rate required to reproduce f and needs to be 
included when comparing with estimates of source emis- 
sivity. 

IV. CONCLUSIONS 

The Lya forest is emerging as a powerful probe of 
structure formation at high redshift. Because the physics 
is relatively well understood we can compare increasingly 
realistic calculations with ever better data to constrain 
our cosmological models. Here we have investigated two 
ways in which the Lya forest constrains cosmology. The 
first, speculative, proposal was to search for baryonic os- 
cillations in the matter power spectrum through their 
imprint in the Lya forest cross-spectra. The second pro- 
posal was to constrain the cmissivity of the sources of the 
UV background using the measured mean optical depth 



of QSO spectra. 

In the first case we found that while in principle the 
signal-to-noise and statistics of the signal were favorable, 
in practice the very low signal levels made this measure- 
ment prone to numerous systematic effects. In the sec- 
ond case we showed that reasonable extrapolations of the 
QSO luminosity function to high-z would allow QSOs 
to contribute a non-negligible fraction of the UV back- 
ground even as early as z — 6. If this turns out to be 
the case, the UV background should have relatively large 
spatial fluctuations which can affect the 2-point statistics 
of the forest at high z. 
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